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ABSTRACT 

Accretion disks around young stars evolve in time with time scales of few million years. We present here a study of 
the accretion properties of a sample of 35 stars in the ~ 3 million year old star-forming region a Ori. Of these, 31 
are objects with evidence of disks, based on their IR excess emission. We use near-IR hydrogen recombination lines 
(Pa7) to measure their mass accretion rate. We find that the accretion rates are significantly lower in a Ori than in 
younger regions, such as p Oph, consistently with viscous disk evolution. The He I 1.083 pm line is detected (either in 
absorption or in emission) in 72% of the stars with disks, providing evidence of accretion-powered activity also in very 
low accretors, where other accretion indicators dissapear. 
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1. Introduction 

, At the time of their birth, circumstellar disks are present 

■ around most (if not all) stars of mass lower than few so- 
I lar masses. The fraction of stars with disks decreases with 
. time, from being close to 100% in the youngest regions to 
' less than few percent after ~ 5 — 7 Myr (e.g., Hernandez et 
[ al. 120071 and references therein). By the time a region is 10 
. Myr old, all the "classical" (i.e., gas-rich, optically thick) 
' disks have disappeared. On similar timescales, all other in- 
, dications of accretion-powered activity, such as accretion 

■ and ejection of matter, disappear as well (e.g., Kenyon et 
I al. 120051 Mohanty et al. 2005; Barrado y Navascues and 
, Martin 2003) . These time scales are consistent to zero or- 

■ der with viscous disk evolution (Hartmann et al. I1998p . 
I although many other processes may play a role in disk dis- 
. sipation (e.g., HoUenbach et al. I2000p . 

Recently, the discussion on how disks evolve has gained 
new momentum from the large surveys of star forming re- 
gions obtained with Spitzer, which provide measurements 
of the fraction of stars with disks based on well-defined, sta- 
tistically significant samples. Particularly interesting is the 
new class of "evolved" disks, of which very few were known 
from ground-based photometry (Skrutskie et al. I1990p . 
They are objects with IR excess emission weaker than typ- 
ical classical T Tauri stars (CTTS); a fraction of these are 
so-called transitional disks, i.e., objects with no near-IR ex- 
cess, but excess emission at longer wavelengths. The frac- 
tion of evolved disks is larger in older regions (Hernandez 
et al. I2007p . The weak emission of evolved disks may be 
the result of grain coagulation and settling (D'Alessio et 
al. HUUBl DuUemond and Dominik HHOS) ■ Transitional ob- 
jects are likely disks where small grains have been cleared 
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from the inner regions, and the discussion on their relevance 
for disk evolution (do all disks disappear from inside out?) 
is open. 

Another indication of the presence of disks is given by 
the accretion-powered phenomena that go under the label 
of accretion activity (UV veiling, line emission and absorp- 
tion from infalling and ouflowing gas, etc.). From them, it 
is possible to measure the mass accretion rate from the disk 
onto the central object, a crucial quantity for constraining 
the physical properties of accretion disks. Recently, Sicilia- 
Aguilar et al. (|2006p measured the mass accretion rate in a 
number of stars of different age from Taurus, p Oph, Cha I, 
TW Hya and Tr 37 and found that Mace decreases steadily 
with time, roughly as t~^-^. Such a trend is consistent with 
the expectations of viscous disk evolution (Hartmann et 

The sample used by Sicilia-Aguilar et al. (|2006p is based 
on individual stellar ages and includes only few stars in 
each region, with the exception of Tr 37, which contributes 
more than half of the total. The mass range of the stars 
with measured accretion rate in Tr 37 (which has a dis- 
tance of ~ 900 pc) is shifted toward higher mass than the 
other, closer regions. We think that it is important to mea- 
sure mass accretion rates in large numbers of stars in dif- 
ferent regions over a large range of average ages. This is 
because there is evidence that, in addition to the time de- 
pendence, other parameters, not identified so far, are lik- 
ley to affect disk evolution, as shown by the MuzeroUe et 
al. ([2003P and Natta et al. (2006) studies of accretion in 
Taurus and Ophiuchus. These authors have shown that, 
firstly, the accretion rate depends on the mass of the central 
object and that, at the same time, there is a large range of 
disk and accretion properties even for stars of similar mass 
and age in the same star forming region. Moreover, there is 
evidence that the disk lifetime is shorter for more massive 
stars (Hillenbrand et al. I1998j Carpenter et al. I2006| Lada 
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We are interested in extending the studies of the accre- 
tion properties of young stars to older regions, to see how 
they change with time and, more specifically, if both the 
correlation between accretion rate and central mass and the 
large spread of accretion rates remain unchanged in time. 
In this paper, we present the first results of a study of the 
accretion properties of CTTS in a Ori. The cr Ori cluster is 
located at a distance of ~ 350 pc and has an age of ~ 3 Myr. 
It contains more than 300 stars, ranging in mass from the 
bright, massive multiple system a Ori itself (spectral type 
09.5) to brown dwarfs. The a Ori region has been exten- 
sively studied over the last few years in the optical. X-ray 
and infrared (e.g., Kenyon et al. 120051 Zapatero-Osorio et 
al.[2002, Franci osini et al. HQDSl Olive ira et al.f2004, 20051 
Caballero 120071 Caballero et al. l2007p . Spitzer has recently 
produced a complete census of the stellar population down 
to the brown dwarf regime (Hernandez et al. I2007p . These 
authors find that the fraction of stars that retain classical 
(i.e., flared, gas-rich) disks (i.e., bona- fide Class II objects) 
vary from about 10% for stars with mass > 2 Mq to about 
35% for T Tauri stars (TTS) and brown dwarfs (BDs). A 
relatively large fraction of TTS, about 15%, have evolved 
disks; 7 of these are candidate transitional disks. The large 
and well characterized sample of Class II and evolved disks 
in cr Ori makes it particularly suited to study the proper- 
ties of accretion and disks for an intermediate age stellar 
population. 

We determine the accretion rate from the luminosity of 
the IR hydrogen recombination lines, Pa/3 and Pa7 in par- 
ticular. The reliability of this method has been discussed by 
M uzeroll e et al. P^M]) . Natta et al. (50D¥) and Calvet et 
al. (|2004p and applied by Natta et al. (2006) to a complete 
sample of IR-selected Class II stars in the young, embed- 
ded star forming region p Oph. Although for the optically 
visible objects in a Ori other methods of measuring the ac- 
cretion rate could (and should) be used, we consider that 
the comparison between different regions is more significant 
if the same method is used in all cases. 

We discuss the properties of our a Ori sample in §2, 
which also presents details of the observations and data 
reduction. The results are presented in §3. The accretion 
rate measurements and their implications are discussed in 
§4. Summary and conclusions follow in §5. 



2. Sample, observations and data Reduction 

2.1. The sample 

Our sample contains 35 TTS with evidence of accretion 
and/or disks. We selected them from two different sources. 
The first is Oliveira et al. (p004l [2006]) . who detected a K-L' 
color excess, very likely due to a circumstellar disk, in about 
25 a Ori TTS. The other source is the high spectral resolu- 
tion spectral survey of Sacco et al. (|2007p . who detected a 
total of 28 stars with broad Ha lines (10%FW>200 km/s), 
indicative of accretion. We imposed a limiting magnitude 
J<14.5, to achieve the required signal-to-noise. Our sam- 
ple includes 18 of the 23 Oliveira et al. sources above the 
J limit, and 22/28 of the broad Ha stars; 5 stars (number 
#5, #6, #8, #14 and #22 in Table[T|) are in common (they 
have both K-L' color excess and strong Ha), for a total of 
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H-[8] 

Fig. 1. [J-H] versus H-[8] color diagram for a Ori cluster 
members with J< 14.4. Empty circles and squares are Class 
Il/transitional disks and Class III objects respectively. Red 
filled circles show the location of our a Ori sample with 
Spitzer observations (32 objects in total). The histograms 
at the bottom show the distribution in H-[8] of all cluster 
members (empty) and of our sample (dashed) , respectively. 
All Class II disks have H-[8] > 1.5. 

35 objects. Sacco et al. p007p confirm the membership of 
the 22 stars in their sample; none is found to be binary. 

After completion of our observations, Hernandez et 
al. (|2007p . pubhshed the results of their Spitzer survey of 
the cr Ori region. Of our 35 objects, 32 lie in the surveyed 
field. Based on the shape of their spectral energy distribu- 
tion (SED) in the near and mid-IR, 26 of them are classified 
as Class II, or classical TTS (CTTS), i.e., objects with gas- 
rich, flared disks extending close to the star; 4 are Class 
HI (i.e., TTS with no evidence of disk) and 2 (#28 and 
#29) are evolved disks, i.e., objects with IR excess weaker 
than typical CTTS. In the following, we will analyze them 
together with the Class II objects and use the definition of 
Class II for the whole disk sample (but see §4.4). The SED 
classification is given in Table [1] Three stars, #1, #13, 
and #33, are not in the fields observed by Hernandez et 
al. f|20D7j) . We classify them as Class II based on the pres- 
ence of a K-L' excess (Oliveira et al. 120041 [2006P : in one case 
(#1) the excess is at the 2.5cr level, and the classification 
uncertain. 

Within our J-magnitude limit (J<14.5), our sample of 
disk objects contains about 30% of the a Ori members with 
evidence of disks from the Spitzer survey. Their distribu- 
tion in color and magnitude is very similar to that of the 
parental population (see Fig. [1]), so that we consider our 
total sample of 31 Class II objects to be representative of 
the whole disk population. 

2.2. Stellar parameters 

Spectroscopically determined spectral types exist for a 
small fraction of cr Ori objects only (e.g., Franciosini et 
al. 2006 and references therein). To ensure homogeneity, 
we determine the effective temperature of each object from 
its visual photometry, adopting the relation between color 
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Table 1. Stellar Parameters and Observed Properties 
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Notes: a, b: Transition disk and evolved disk, respectively, in Hernandez et al. (|2007p . c: tentative {2a) detection, d: Type IVr 
P-Cygni profile: EW(em) = 1.2±0.3, EW(abs)=-0.9±0.3. e: tentative type IVr P-Cygni profile (0.6, -0.6). 



indexes, T^ff and spectral types of Bessell ()199ip for M 
dwarfs stars, and from Bessell (|1979p and Bessell and 
Brett (p^988) for stars of earlier spectral type. We assume 
negligible extinction in all bands (Brown et al. I1994t Bejar 
et al. T999) , and compute Teff independently from the 
three color indexes (V-R), (V-I), (R-I). The optical pho- 
tometry is from Sherry et al. (|2004p . Zapatero-Osorio et 
al. ([2 ^1), K enyon et al. p005)) . Bejar et al. pOOT)) and 
Wolk p996p . The differences between the different deter- 
minations of Teff are generally less than ± 150 K, and 
often much smaller. The adopted values of Teff and spec- 
tral types are given in Table [TJ The comparison with the 
spectroscopic determinations of spectral types of Zapatero- 
Osorio et al. (|2002p for the 8 objects in common shows 
differences of 1/2 to 1 subclass. 

Stellar luminosities have been computed from the I mag- 
nitude, using the bolometric correction appropriate for the 
spectral type (Kenyon and Hartmann '19950 , no extinction 
and a distance D=350 pc. The location of the stars on the 
HR diagram is shown in Fig.[2l together with the evolution- 
ary tracks of D'Antona and Mazzitelh (DM97: [1997^ . The 



distribution of the stars is consistent with a typical age of 
~ 3 Myr (Oliveira et al. l2004p . with the exception of #9, a 
Class III object with color indexes consistent with an early 
K spectral type but very low luminosity. 

The stellar masses range from ^ 0.5 to ^ 0.1 M0. Note 
that the exact values of depend on the adopted evolu- 
tionary tracks, especially around 0.5-1 Mq. We use DM97 
for homogeneity with the analysis of Natta et al. (|2006p of 
the mass accretion rates in p Oph. The lower mass limit 
is determined by our selection criterion that J<14.5 mag; 
the lack of more massive stars reflects the strong decrease 
of the fraction of Class II objects for increasing J mag- 
nitude observed by Hernandez et al. ()2007 |). This trend, 
seen in other star forming regions (e.g., Hillenbrand et 
al. 1998; Carpenter et al.|2006; Lada et al. [20061 Hernandez 
et al. l2007j Dahm and Hillenbrand I2007p . is particularly 
strong in a Ori, where only 4 (7%) stars brighter than J=ll 
mag are classified as Class II. One of them is included in 
our sample (#23), and has an estimated spectral type K8 
and mass ~ 0.5 Mq. The fact that no higher-mass member 
has an associated Class II disk is confirmed by Caballero 
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Fig. 2. Location of the a Ori targets on the HR diagram. 
Evolutionary tracks from D'Antona and MazziteUi (1997|, 
labeUed with the steUar mass. Isochrones for 0.5, 1 and 5 
Myr are shown by the dashed hnes. 

(|2007p : of the 18 very bright stars (spectral type GO and 
earlier) in the a Ori cluster studied by him and included 
in the list of cluster members by Hernandez et al. (|2007p . 
3 have transitional or debris disks, all the others are Class 
III. 

2.3. Observations 

Near-infrared spectra in the J band were obtained for all our 
targets using SOFI near-infrared camera and spectrograph 
at the ESO-NTT telescope. The observations of a Ori were 
carried out in Visitor Mode on December 1-4, 2006. The 
0.6 arcsec slit and the Blue low resolution grism have been 
used, resulting in a spectral resolution A/AA ^ 1000 and 
a spectral coverage from ~0.95 to ~1.64/im. Integration 
time was 40 minutes for all the targets; standard calibra- 
tions (flats and lamps) and telluric standards spectra were 
obtained for each observation. Wavelength calibration was 
performed using the lamp observations. 

The mean seeing during the three nights was 0.7 arc- 
sec, with many hours below 0.5 arcsec. All target objects 
and telluric standards have been observed at airmass values 
z<2. 

With the goal of observing as many objects as possible, 
the targets have been observed in pairs, choosing stars with 
angular separations which allow a 20 arcsec nodding and 
a 6 arcsec jitter with the 4.71 arcmin slit, and orienting 
the slit at appropriate position angles. In total, we have 
observed in this way 24 targets, while 11 have been observed 
individually. 

The spectra have been reduced using standard proce- 
dures in IRAF. We are interested in this paper in the three 
lines indicative of accretion powered activity which fall in 
the J band, namely the Hel line at 1.083 ^m and the two 
hydrogen recombination lines Pa7 at 1.094 //m and Pa/? at 
1.282 /zm; the portions of the spectra which contain these 
lines are shown in Fig. [3] and [J] Their equivalent widths 
are given in Table [1] (positive values for emission lines). 
Upper limits to Pa/? and Pa7 assume that the lines are in 



emission; they have been estimated from the peak-to-peak 
noise in the spectral region of the lines, under the assump- 
tion that the lines are not resolved. Since the He line can 
be either in emission or in absorption or both, we do not 
explicitely give upper limits in Table (TJ they are similar to 
those of the adjacent Pa7 line. 



3. Results 

The presence of the hydrogen lines in emission and of the 
He I line, often with P-Cygni profiles clearly seen in high 
resolution spectra, (Edwards et al. 2006) is evidence of ac- 
cretion related activity. In fact, inspection of Table [T] shows 
that none of the three lines is detected in Class HI stars, 
with the possible exception of #15, where we have a ten- 
tative detection of He I in emission. Of the 31 Class II, we 
detect at least one of the three lines in 26 cases. The hydro- 
gen lines are detected in emission in 12 Class II stars, but 
only in 5 we do detect both of them; in 7 cases, only Pa7 
is clearly seen. The hydrogen line detection is confirmed 
in 10/12 cases by the presence of the He I line, either in 
absorption or in emission. 

Our spectral resolution 300 km/s) is too low to re- 
solve the lines. However, some lines appear to be marginally 
resolved. It is possible that weak lines look broad due to 
noise spikes in the adjacent continuum, as we think is the 
case of Pa7 in object #3, which has an unrealistic FWHM 
of 950 km/s; we did not try to correct for this effect the 
equivalent width in Table [U which is hence very likely over- 
estimated by a factor of ^2. In two stars (#19 and #21), 
the He I line has an inverse P-Cygni profile, with redshifted 
absorption; the line is strong in #19, but barely detected 
in #21. Table [T] gives the total equivalent width of the line, 
while the values of the absorption and emission components 
are given in the notes to the table. 

We have been surprised by the number of objects where 
Pa7 is detected while Pa/3 is not. The ratio of the Pa7 to 
Pa/3 fluxes has been recently measured for a group of CTTS 
in Taurus by Bary and Matt (2007), who find a tight cor- 
relation with a ratio 0.86±0.10. Magnetospheric accretion 
models (Muzerolle et al. l200ip predict ratios ^ 0.7 — 0.9 for 
MaccZ10~^ Mq/jt, increasing to 1.2-1.4 at Macc^ 10"^ 
Mq/jt (J. Muzerolle, private communication). Considering 
that for TTS the ratio of the continuum at the two line 
wavelengths is ~ 1.2 — 1.4, we expect in any case compa- 
rable equivalent widths, and, in fact, this is the case in the 
five a Ori objects where both lines are detected. All stars 
with Pa7 but no Pa/3 have relatively weak Pa7; however, 
in several cases we should have detected Pa/3 if the equiva- 
lent widths were similar and the lines unresolved. We have 
checked the observational and data reduction procedures, 
and found no obvious explanation for the absence of Pa/5. 
It is possible that the upper limits to Pa/3, derived assum- 
ing that the line is not resolved, are in fact underestimated, 
and, indeed, some objects with very low limits to the ra- 
tio of the Pa/3/Pa7 equivalent width (e.g., #3, #22, #26) 
have broad Pa7 (see Fig. [3]). Also, in some cases Pa/3 may 
have a P-Cygni profile with deeper absorption than Pa7, 
and this could reduce the ratio of the unresolved equiva- 
lent widths. This kind of profile, however, is quite rare in 
Tauru s TTS (e.g., Edwards et al. 2006, Folha & Emerson 
I200ip . Note that the two lines Pa7 and Pa/3 are observed 
simultaneously, so that variability cannot affect their ratio. 
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Fig. 3. Normalized observed spectra in the Hel and Pa-y region. The dotted hnes show the wavelengths of the two lines. 



The Hel 1.083/im line is seen in 23/31 Class II stars, in 
absorption (8 cases), in emission (14 cases), or with a P- 
Cygni profile (2 cases). This variety of profiles is known to 
occurr in TTS (Edwards et al. I2006p . and will be discussed 
in S 4.2. 



4. Discussion 

4.1. Accretion luminosities and mass accretion rates 

We use the hydrogen recombination lines to measure the ac- 
cretion luminosity Lace, following the procedure described 
in detail by Natta et al. (pOng]) . 



Firstly, we compute line fluxes from the equivalent 
widths and the flux of the continuum near each line, cal- 
ibrated with the 2MASS value of the J magnitude and 
Aj=0. We apply a correction of a factor 1.3 to the con- 
tinuum at 1.09 yum, to account for the typical slope of the 
J-band continuum in late K - M stars. 

Line luminosities are computed from the line fluxes for 
the adopted distance D=350 pc. We use them to derive 
a measurement of the accretion luminosity Lace, using the 
empirical correlation between the luminosity of the near- 
IR hydrogen recoombination lines and L^cc first noticed 
by MuzeroUe et al. (|1998p and used to study the accretion 
properties of TTS in Ophiuchus by Natta et al. pOOe]) . 
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Fig. 4. Normalized observed spectra in the Pa/3 region. The wavelengtth of Pa/3 is shown by the dotted hne. 



The correlation has been established quantitatively for the 
two lines Pa/? and Br7 in a sample of TTS in Taurus for 
which reliable measurements of the accretion luminosity 
from veiling were available (MuzeroUe et al. 119981 Calvet 
et al. 120041 , and extended in the case of Pa/3 to the brown 
dwarf regime by Natta et al. (|2004p . We adopt here a re- 
lation between Lqcc and Pa7 luminosity obtained from the 
Natta et al. (|2004p correlation between Lace and the Pa/3 
luminosity and the ratio of Pa7 to Pa/3 fluxes of 0.86 ±0.1 
(Bary and Matt IM?7|) : 

Log^ = 1.36 X Log^^^^ + 4.1 (1) 



Eq.(l) does not depend on any assumption on the line 
origin, but only on the known empirical relation between 
Lqcc and the Pa/3 luminosity and the adopted value of the 
Pa7/Pa/3 ratio. As discussed at the end of §3, in several of 
our a Ori objects the observed ratio Pa7/Pa/3 is higher than 
the Bary and Matt one and the values (or upper limit) of 
Lqcc computed from L(Pa/3) are correspondingly lower (see 
Fig. [5]). We decided to use in the following Lace derived from 
Pa7 to maximize the number of actual measurements; since 
Lacc(Pa/3) < Lacc(Pa7), deriving Lace from L(Pa/3) would 
just strengthen our conclusions. 
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measurement errors on the line equivalent widths, there are 
continuum calibration uncertainties, since our spectra have 
not been photometrically calibrated. Variability in the lines 
is a well known phenomenon in TTS, and snap-shot esti- 
mates of their strength can differ from the long-term aver- 
age values by large factors. Moreover, the relation between 
L(Pa7) and Lqcc has an intrinsic, non-negligible scatter, 
and the determination of the stellar mass/radius depends 
on a number of assumptions, among them the adopted evo- 
lutionary tracks. The final uncertainties on Mace can easily 
be of ±0.5 dex, as discussed, e.g., by Calvet et al. (|2004p 
and Natta et al. (2006). In spite of all that, snap-shot ac- 
cretion values of a large number of stars in any given star 
formation region provide very significant information on the 
typical accretion properties of the region and on its evolu- 
tionary stage. 

The results for a Ori are summarized in Fig. [5] and [71 
top panel. Fig. [5] plots the accretion luminosity as a func- 
tion of the stellar luminosity for all Class II objects. It shows 
clearly that very few stars have L^cc larger than 0.1 L.^, and 
that most of them have values well below this limit. This 
is different from the results obtained in younger star form- 
ing regions, such as Taurus and Ophiuchus (see summary 
in Natta et al. I2006p . where a large fraction of stars has 
IU> 0.1. 

Fig. [7] shows the mass accretion rate as function of the 
stellar mass for the same 31 objects. The values of Mace 
are all smaller than ~ 10~* Mo/yr, and 23/31 stars have 
accretion rates lower than 10~^ Mo/yr. When compared to 
stars with similar mass in Taurus and Ophiuchus, these are 
very low values, as we will discuss in § 4.5. 
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minosity for the Class II objects. Dots are actual measure- 
ments, arrow 3ct upper limits. The dotted lines show the 
location of the Lqcc/L* = 1-0, 0.1 and 0.01 ratios, respec- 
tively. 



The mass accretion rate Mace is then been computed 

from Lace- 

where M, and R» are the stellar mass and radius, respec- 
tively. 

The uncertainties on Lqcc and Mace are difficult to es- 
timate, but they are certainly not small. In addition to the 



4.2. He line 



Edwards et al. (|2006p have recently published a study of the 
Pa7 and He I profiles in a sample of 39 Taurus TTS; they 
found that the strength of both lines correlate with the J- 
band veiling, i.e., with the accretion luminosity. Although 
with a large spread, objects with low veiling have on av- 
erage not only weaker but also narrower Pa7. He I lines 
have very often P-Cygni profiles, with blue-shifted absorp- 
tion caused by outflowing gas and red-shifted absorption 
due to infalling material. On average, strong accretors have 
the He I line in emission; red-shifted and blue-shifted ab- 
sorption features are present in high and low accretors, but, 
if unresolved, the He I equivalent width would turn from 
strong emission in high accretors to net absorption in low 
accretors. In these objects, the blue-shifted absorption fea- 
tures due to outflowing matter, very likely a non-coUimated 
wind, dominate the profile. 

Fig. [5] plots the equivalent width of the He I line as 
function of the accretion luminosity L^cc for our sample. 
If we consider objects with measured hace, we find a good 
correlation between the two quantities, with He I in net ab- 
sorption in low accretors and in emission in high accretors, 
as in the Edwards et al. (2006) sample. Note that Edwards 
et al. use in their analysis the "activity index" of the He I 
line, i.e., the sum of the absolute values of the emission and 
absorption equivalent widths, while our unresolved profiles 
give the net equivalent width only. The trend of decreasing 
emission, increasing net absorption for decreasing Lace is 
weaker if we include objects with Lace upper limits; as in 
Taurus, there are some objects with strong He I emission 
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but low Lcicc(e.g., #27), which it would be interesting to 
investigate further. 

The Taurus sample includes CTTS with mass accretion 
rates Mace between '--^ 5 x 10^^ and 5 x 10^^" Mo/yr, with 
a large fraction (29/38) of objects accreting at rates higher 
than the highest values of the a Ori sample. Edwards et 
al. (|2006p divide their sample of CTTS into four groups ac- 
cording to the 1/im veiling. The median values of Log Mace 
of the first and second group (high and medium veiling) 
are -5.8 and -7 Mo/yr, while the third and fourth (low 
1/im veiling, narrow and broad Pa7, respectively) have me- 
dian -8.2 Mo/yr. The low veiling, broad Pa7 group has, 
on average. He I lines with net negative equivalent widths. 
The a Ori sample overlaps with the low veiling objects and 
extends it to even lower values of Mace- In fact, there are 
no objects with the very strong He I emission lines of some 
Group I Taurus stars, while the fraction of objects with net 
He I absorption is higher than in the Edwards et al. (|2006p 
group IV. 

An interesting aspect of these results is that absorp- 
tion He I lines can be detected even in objects with very 
low Maccj where direct accretion signatures (e.g., hydrogen 
line emission, veiling) can be hard to measure. Accretion- 
powered winds may be the longest lasting tracers of disk 
activity. 

4.3. Evolved Disks 

The a Ori cluster contains a relatively high fraction of 
evolved disks, i.e., objects with IR excess emission lower 
than the median of CTTS. Some of them are classified by 
Hernandez et al. (|2007p as transitional disks, i.e., objects 
with a very low excess in the near-IR but a normal excess at 
longer wavelengths. Our sample includes two such objects. 
One (#28, S897) is classified as a transitional disk; it has no 
Pa/3 or Pa7 lines above the detection limit (the correspond- 
ing mass accretion rate is < 10~^ Mo/yr), but a tentative 
detection of He I, and a rather broad Ha (10%FW=504 
km/s; Sacco et al. IMTTI) . The other object (#29, S908) is 
classified as an evolved disk, but shows both Pa/3 and Pa7 
in emission and the He I line in absorption. Also for this 
object, Sacco et al. (2007) measure a broad Ha, with a 
10%FW=401 km/s. Its accretion rate is similar to that of 
the Class II a Ori stars. 

It is possible that in these objects (#29 in particular) 
the low IR excess is due to the disk orientation on the plane 
of the sky. However, some transitional disks are known to 
have relatively large Mace (see summary in Chiang and 
Murray-Clay pOOTj) and references therein), and the dis- 
cussion on how long accretion can be sustained once disks 
start to evolve is still open. The relatively large sample 
of evolved and transitional disks in a Ori is well suited 
to study the accretion properties of transition and evolved 
disks. However, the sample discussed in this paper includes 
too few such objects to provide any answer. 

4.4. Accretion Evolution 

The relevance of the a Ori results for the understanding of 
the accretion evolution can be better understood in com- 
parison with other regions of star formation, both younger 
and older. 




Log M^/Mq 

Fig. 7. Mass accretion rate as function of M* for a Ori 
(top) and p Oph (bottom; from Natta et al. l2006p . Dots are 
actual measurements, arrows 3cr upper limits. The dotted 
line (same in both panels) shows the log Mace = —7.3 -t- 
1.8 X logM^, locus, that we take as a fiducial line to define 
"high accretors" (see text). 



The comparison with the p Oph star forming region is 
particularly interesting. Mass accretion rates in p Oph have 
been derived by Natta et al. (|2006p for an IR-selected sam- 
ple of Class II complete to a limiting mass of about 0.03 
Mo, using the luminosity of near-IR hydrogen recombina- 
tion lines as a proxy of the accretion luminosity. Because 
the selection criteria and methods are similar, the compar- 
ison between these two regions is more meaningful than in 
other cases. 

The major difference between the two regions is the lack 
in a Ori of objects with high accretion rates. If we restrict 
the comparison to the mass range 0.1-0.5 Mo, where the 
two samples overlap, and take as a reference the line shown 
in Fig. [71 which is the best fit to the p Oph sample of 
measured Mace over the whole mass interval, we find that 
only 1 star in cr Ori lies above this line (3%), while the 
fraction is ~30% in p Oph. 

Because the a Ori sample is small (much smaller than 
the p Oph one), we performed a simplified Monte Carlo 
simulation to assess the statistical significance of this dif- 
ference. The goal of our simulation is to evaluate whether 
the number of high accretors we observe in the a Ori sam- 
ple is consistent with being drawn from a distribution of 
rates as derived from the p Oph Class II objects observa- 
tions. We take the p Oph stars in the interval 0.1-1 Mo 
(the sample 0.1-0.5 Mq would have too few sources to de- 
rive a reliable Mace distribution) and use the relationship 
between accretion rates and stellar mass derived by Natta 
et al. (j2006p to scale all measurements to a fiducial stel- 
lar mass of 1 Mq; we then compute the fraction of high 
accretors as the ratio of objects with measurements above 
the fiducial line shown in Fig. [7] to the total (measurements 
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plus upper limits). Within the chosen mass range, the frac- 
tions of high accretors in p Oph and a Ori are 29% and 3% 
respectively. 

Using a Monte Carlo method, we extract a large number 
of samples of 31 accretion rates. For each sample we ran- 
domly remove a number of objects based on the fraction 
of upper limits in p Oph (~52%), while for the remaining 
objects we simulate measurements using the distribution of 
measured rates. Figure [9] summarizes the results: the his- 
togram plots the distribution of the fraction of high accre- 



tors in the simulations, the dashed line shows the value for 
the p Oph sample while the dashed-dotted line shows the 
value observed in our a Ori sample. The value observed in 
a Ori is lower than the 0.1 percentile of the Monte Carlo 
simulations, suggesting that it is very unlikely that the two 
samples are two different realizations of the same popula- 
tion. 

The lack of strong accretors in a Ori is very likely an 
age effect, as in viscous disks the mass accretion rate is 
expected to decrease with time (roughly as t~^'^; Hartmann 
et al. ll998p . On average, this means approximately a factor 
10 from the young (less than 1 Myr) p Oph population to 
the older (~ 3 Myr) a Ori one, consistent with the results 
shown in Fig. [71 One way of looking at this is to compute 
the fraction of high accretors in a Ori not with respect 
to the "younger" p Oph Macc-^* fiducial fine, but with 
respect to an "aged" version, where Mace has decreased by 
the same factor 10 for each M*. Then, the fraction of a Ori 
high accretors will be of ^30%, identical to the fraction of 
high accretors in p Oph. 

Fig. [TU] shows a comparison between mass accretion 
rates in a Ori and in the ^--^4 Myr old star forming re- 
gion Tr 37. This region, located at a distance of ^ 900 pc, 
has been extensively studied by Sicilia-Aguilar et al. (20011 
120051 l2006l I2006p . who have classified disk properties from 
Spitzer IRAC and MIPS data and measured Mace from U 
band excess and Ha profiles. The methods of deriving Mace 
are different in the two regions, but probably any system- 
atic effect is within the uncertainties, since both correla- 
tions between the U-band excess and Lace and between the 
hydrogen IR line luminosity and hace have been empiri- 
cally calibrated using the same sample of TTS in Taurus 
with good measurements of hace from veiling (e.g., Calvet 
and GuUbring fTMSl MuzeroUe et al. fTOMl) . 
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Fig. [TO] plots all the stars with disks (bone- fide Class II 
and transition disks) with measured stellar mass and Mace, 
as reported by Siciliar-Aguilar et al. (2006,) . We have not 
included any of the very young (~ 1 Myr) Class II objects, 
probably an episode of triggered star formation, found in 
the TR 37 "globule". Because of its larger distance, the 
Tr 37 sample covers a mass range between ^ 0.4 and 1.8 
M0, with one star of 2.5 (whose Mace is very uncertain); 
there is practically no overlap with the a Ori lower-mass 
sample. However, the accretion properties of the two regions 
look very similar. The trend of higher Mace for higher M, 
is confirmed, with good merging between the two samples. 
More importantly, both a Ori and Tr 37 lack high accretors 
to a similar degree. If we define, as before, high accretors 
to be those stars with values of Mace above the fiducial line 
s hown in Fig. [3 the fraction of high accretors in Tr 37 
is 6%, similar, within the uncartinties, to the 3% found in 
a Ori, and certainly much lower than the 30% fraction of 
high accretors in p Oph. 

The continuity in the Macc-M, behaviour between the 
two regions seems to extends to stars more massive than 
about 1-1.2 Mq. In other star-forming regions, as discussed 
in §1, stars more massive than ^ 1 Mq loose evidence of 
disks and accretion much faster than less massive stars, and 
one could expect to see a drop in the fraction of Class II 
stars among the more massive members of Tr 37. However, 
this does not seem to be the case, because Tr 37 has a 
relatively high fraction of Class II (8 with respect to 5 Class 
HI among stars more massive than 1.2 Mq), even when the 
globule population is not taken into account. Of these, 6 
have evidence of accretion. The case of Tr 37 shows that 
the evolution of disks and accretion as function of the stellar 
mass is still an open problem, that needs to be addressed 
in more detail in the future. 

5. Summary and conclusions 

This paper presents J-band spectra of 35 stars in the star- 
forming region a Ori. We concentrate our attention on the 
three lines that are indicative of accretion-related activity, 
namely the two hydrogen recombination line Pa/3 and Pa7 
and the He I line at 1.083 ^m. 

The stars range in mass from ^ 0.1 to ^ 0.5 Mq. Four 
have no evidence of disks (Class HI), according to the classi- 
fication of Hernandez et al. (|2007p based on Spitzer near and 
mid-IR photometry. None of the three lines is detected in 
their spectra, with the possible exception of one star which 
has a tentative detection of He I emission. The other 31 
objects have associated disks (Class II and evolved disks), 
and we detect Pa7 in 12 of them and Pa/3 in 5. The He I 
line is seen in 23 disk objects, either in emission or in ab- 
sorption or both. Even at the low spectral resolution of our 
spectra (~ 300 km/s), in 2 cases the He I has a P-Cygni 
profile with redshifted absorption. 

We derive accretion luminosities and mass accretion 
rates from the Pa7 luminosity; the results for a Ori indi- 
cate values much lower than in stars of similar mass in the 
younger regions Taurus and p Oph. In particular, a statis- 
tical analysis shows that the p Oph and a Ori populations 
cannot be drawn from the same parent population. TTS in 
a Ori are statistically identical to those in p Oph only if 
we take into account their age, i.e., that any p Oph object 
will have, at a Ori age, a value of Mace lower by a factor 



of about 10 {Maccoc t^^'^), as expected if viscosity controls 
disk evolution. 

The Hel is detected (in net absorption or emission) 
in more stars than either of the hydrogen lines. Its be- 
haviour agrees qualitatively with what found by Edwards 
et al. (|2006p in a sample of TTS in Taurus; namely, strong 
Hel emission characterizes on average high accretors, while 
net absorption, probably due to a wind, is typical of low ac- 
cretors. We extend this result to objects of very low Mace 
and suggest that the He I line may be used to detect weak 
accretors when other methods become unfeasible. 

The a Ori region, with its well studied disk population 
and wealth of data at all available wavelengths, large num- 
ber of member stars extending over a broad mass range and 
moderate distance, is particularly well suited for a number 
of follow-up studies. Among them, we want to mention that 
it will be relatively easy to obtain near-IR spectra of the 
complete sample of Class II objects, down to the brown 
dwarf regime, to check, e.g., if the accretion rate evolves 
differently in brown dwarf disks. The comparison of a Ori 
with the region Tr 37, which has similar age and a com- 
parable fraction of high accretors, but a surprisingly high 
fraction of disks among stars more massive than about 1 
Mq , shows that the dependence of disk and accretion evo- 
lution with the mass of the central object is still an open 
and intriguing problem. 

Our results on the two evolved/transition disks in our 
sample confirm that accretion is still going on in some of 
them and suggest that near-IR spectroscopy and high res- 
olution spectra of the hydrogen and He I lines may help to 
clarify the last stages of the history of accretion disks. 
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